Measuring heavy neutrino couplings at the LHC 
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The existence of heavy neutrinos mediating neutrino masses via a type-I seesaw can be directly 
probed at the LHC, or indirectly in rare lepton flavor-violating processes. The synergy between 
these two approaches requires a direct measurement of the heavy neutrino couplings. We discuss a 
strategy to perform such measurements at the LHC in the context of left-right symmetric models, 
which is also applicable to other models implementing a type-I or type-Ill seesaw. We demonstrate 
that the ambiguities in the determination of the heavy neutrino mixing parameters can be resolved 
by performing an exclusive analysis of dilepton final states, discriminated by flavor and missing 
energy. 



I. INTRODUCTION 

The smallness of light neutrino masses, compared to 
the other fermions, suggests that new physics beyond the 
standard model (SM) is involved in their generation. The 
simplest and most popular explanation for neutrino mass 
suppression is provided by the seesaw mechanism 
which introduces heavy right-handed neutrinos. These 
new states need not to be at the TeV scale, nor they 
are necessarily at the direct reach of the Large Hadron 
Collider (LHC). Nevertheless, this is an interesting pos- 
sibility which has motivated much effort both from the 
theoretical an d experimental [ljj [2(| sides. 

It is well known that the heavy neutrinos may en- 
hance the rates of lepton flavor-violating processes such 
as n — > e7, Z — » efi, /i — > eee and fj, — e conversion 
(see [21], [22j for reviews). In this sense, precision measure- 
ments in the intensity frontier constitute a complemen- 
tary window to testing the seesaw mechanism. Notably, 
the interplay between the high-energy and high-intensity 
probes of the seesaw requires the direct measurement of 
the heavy neutrino couplings at the LHC. In this work 
we investigate how these couplings could be measured in 
case a heavy neutrino signal is observed. 

As a benchmark scenario, we consider heavy neutrino 
production in the context of left-right symmetric mod- 
els for which the LHC reach is very good 0-01 ■ In 
these scenarios, heavy neutrinos can be produced from 
the decay of an SU(2).r gauge boson Wr, 



pp — > Wr — > IN , 



(1) 



and they subsequently decay through an off-shell Wr as 

N^£'jj,£'tb, (2) 

giving a It' jj signal, where t, £' = e, /i, r (with the same 
or opposite charges, due to the Majorana nature of N) 
and j are light jets. Until now, almost all studies 
dedicated to the study of heavy neutrino signals in left- 
right models have been concentrated on the simplest case 
where i,£' = e,fj,, ignoring the production of r leptons. 
However, current experimental searches are sensitive to 



the r's [6] through the detection of secondary leptons 
produced in the leptonic decays r — > e/fiuv. We stress 
that, from the theoretical point of view, there is no com- 
pelling argument in favor of zero tN mixing. Actually, 
bearing in mind the lepton mixing pattern revealed by 
neutrino oscillation experiments, one could easily find a 
motivation against that assumption. In view of this, an 
analysis including the r-induced dilepton final states is 
most welcome. 

In this paper we will assume that a ii'jj heavy neu- 
trino signal is seen at the LHC. After determining the 
Wr and heavy neutrino masses, which is quite straight- 
forward from the observed resonance peaks, the next step 
would be to measure the couplings of the heavy neutrino 
to the charged leptons and the Wr boson. This interac- 
tion is parametrized by the Lagrangian 

C = -^V&Ih^NrW^ + H.c. , (3) 

where gR and V/^ are the SU(2)# gauge coupling and 
tN mixing parameters, respectively. The latter obey the 
general normalization condition 
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WR 
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(4) 

with k constrained to be very close to unity. For simplic- 
ity, in this work we consider that only one heavy neutrino 
N is lighter than the Wr boson. In case there is more 
than one, the method presented here can be extended 
by applying suitable kinematical cuts to single out each 
heavy neutrino contribution. 



II. MEASURING HEAVY NEUTRINO 
MIXINGS IN DILEPTON EVENTS 

Let us now discuss in some detail how the mixings can 
be extracted from event counting in dilepton final states, 
without assuming unitarity of the leptonic mixings, i.e. 
k ^ 1. The cross section for the process in Eq. ([1) with 
a given lepton £ is 



a{£N) = <t Wr 



^ 2 + R 



(5) 
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where ctw R is the cross section for Wr production and 

Thad 



R 



r 



Lep 



(G) 



is the ratio between the hadronic and leptonic Wr decay 
widths for n = 1 (this quantity is independent of lepton 
mixings). The total cross section for N production is 
then 



Otot — &Wr 



K 2 + R 



(7) 



which, together with the branching ratio for the decay in 
Eq. ©, 



Br(A -> £'X) = 



\V R I 2 

I v l'N\ 
K 2 



(8) 



leads to the cross section for the production of charged 
leptons £, £', 



(J tot it/_R |2|t/-R |2 

aw - — \Vw\ \ v i'N\ ■ 

K 

By introducing the normalized mixings 



yR 
R _ V £N 
V IN — 



we find the expected result 

(TM = 'tot|w^ v | 2 |'i;^ A r| 2 , 



(9) 



(10) 



(11) 



i.e., the total heavy neutrino cross section is distributed 
among the different dilepton flavors according to the 
squares of the normalized mixings. The above equation 
also reflects the fact that one cannot simultaneously mea- 
sure the total cross section and test the unitarity of the 
mixings from au' alone. However, if the gauge coupling 
gn could be determined from other processes (for exam- 
ple, in the decay channel Wr — > tb) and a definite predic- 
tion for aw R and R could be made, then the unitarity of 
leptonic mixing would be testable by inverting Eq. (J7J, 



1 — Otot/CWji 



1/2 



(12) 



We now address the question on how the individual 
cross sections aw (and, thus, the mixing parameters) 
can be extracted from the data. We restrict ourselves to 
final-state signatures with two electrons or muons, which 
can result from prompt production £, £' = e, fi or from 
the leptonic decay of £, £' — r. Notice that with the tag- 
ging of hadronic r decays the final state of one electron 
or muon plus a r-jet could be used as well, improving the 
results presented here. In a previous work Q, it has been 
emphasized that the key for the sensitivity to £, £' = r in 
the dilepton final state is to drop any requirement on the 
missing energy of the events, E™ lss , which is small (and 
caused by instrumental effects) for £,£' = e, /i but large 



when there is one or more r leptons. However, this is 
not enough to obtain precise measurements of the mix- 
ings and, therefore, a separation of the dilepton samples 
according to E^ lss is necessary. This fact, which will be 
confirmed later by our numerical results, can be easily 
understood with the following simple semi-quantitative 
argument. Let us assume the same efficiency e for elec- 
tron and muon detection. At high transverse momenta, 
as it corresponds to the production of a heavy Wr, the 
efficiency for secondary leptons originated from r decays 
is similar but these are penalized by the branching 
ratio Br(r — > evv) ~ Br(r — > [ivv) ~ 0.17 = rj. Then, 
for an integrated luminosity L the number of events in 
each dilepton channel is 



N ee = 
= 



Nu 



a tot Ls (kV + 2r,\v? N f\v? N f + r, 2 |^| 4 ) , 
2a tot Le [\vf N \ 2 \v^ N \ 2 
+v(\v? N \ 2 + \v2 N \ 2 )\v? N \ 2 + V 2 \v? N \ i ] , 
^ot^(|< JV | 4 + 2r,|^| 2 K^| 2 +ry 2 |^ JV | 4 ) . 



(13) 

Although there are three independent parameters (for 
example, v^ N , v^ N and ctot) and three measurements, 
it is clear that very distinct parameter space points can 
give roughly the same numbers of events in the three 
channels. 1 In particular, along the line v^ N = v^ N one 
has N ee = JV W = l/2N efl , and the system of equations 
is underconstrained. For instance, the cases of maximal 
eN and fiN mixing with a tot = °o 



,R — „.R 



1 



v eN = V ^N = -j= , V ?N = °> Vtot = 00 , (14) 



and the one with r mixing only and Utot = co/v 2 ? 
•fit 



VeN = v nN=0, V? N = l, CTtot = Volrf , (15) 



lead to the same number of events. The obvious differ- 
ence between these two cases is that, in the latter, the 
events have large missing energy due to the neutrinos 
produced in t decays. Then, the solution to the degen- 
eracy is to separate the dilepton samples in 'low £^ lss ' 
and 'high £?™ ISS ' events, according to some convenient 
threshold. This distinction doubles the number of mea- 
surements and definitely resolves the aforementioned am- 
biguities, which are otherwise present when event count- 
ing is performed on the i?™ lss -inclusive samples. 



III. NUMERICAL ANALYSIS 

In order to test the above method in the extraction 
of heavy neutrino mixings, we rely on a fast simulation 



1 This degeneracy is not found in 0| , because V^ N = is implicitly 
assumed. 
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analysis of heavy neutrino signals and corresponding SM 
backgrounds at the LHC, with a center of mass energy 
of 14 TeV. A luminosity of 30 fb^ 1 is assumed. Signal 
events are generated with Triada [l7[ and backgrounds 
with Alpgen [24 |. The hadronization is performed by 
Pythia [Hi, and AcerDet (H] is used for a fast sim- 
ulation of a generic LHC detector. For the suppression 
of the SM backgrounds we apply an event selection sim- 
ilar to the one used by the ATLAS Collaboration for 7 
TeV searches [2(| , but tightening some of the kinematical 
cuts. Namely, we require: 

• Exactly two leptons with transverse momentum 
pr > 50 GeV, and pseudo-rapidity < 2.47 for 
electrons, \rj\ < 2.4 for muons. Electrons in the 
range 1.37 < |?7| < 1.52 are excluded. 

• At least one jet with pr > 20 GeV and < 2.8. 

• Dilepton invariant mass to« > 110 GeV. 

• The sum of transverse momenta of the two leptons 
and highest pr jets (including up to two jets) must 
be larger than 700 GeV. 

• The Wr reconstructed mass must be larger than 
700 GeV. For events with at least two jets, this 
mass is defined as the invariant mass of the two 
leptons and the two leading jets, M^ Q R = rriujj- 
For events with only one jet, it includes only this 
jet, Mffi R = m Uj . 

We split our dilepton final state events in twelve dis- 
joint sets. First, those with same-sign leptons are sepa- 
rated from the ones with opposite-sign. This turns out 
to be convenient since the former have a much lower 
SM background. Second, each one of those two sets 
is divided in two subsets with E'™ ss < 200 GeV and 
£™ 1SS > 200 GeV, corresponding to 'low' and 'high' 
E™ lss . Finally, the events are classified according to the 
lepton flavour. The SM background Bx for each final 
state generically labelled as X is given in Table U for a 
luminosity L = 30 fb~ . It includes tt, single top, W/Z 
+ jets, W/Zbb + jets, W/Zcc + jets, W/Ztt + jets, 
diboson and triboson production. Further details about 
the background generation can be found in (l7j . 

TABLE I: SM background Bx for the twelve dilepton final 
states X and for a luminosity of 30 fb -1 



Same sij 


;n 


Opposite sign 


Low Eg lss Hi; 


;h Eg lss 


Low E^ iss High E!? iss 


ee efi fifi ee 


efi fifi 


ee efi fifi ee efi fifi 


11 15 8 3 


3 2 


1180 702 1316 91 203 227 



For the heavy neutrino signal we assume g^ — g^ and 
consider two mass points: (1) Myy R = 2.5 TeV, = 1 
TeV, for which o Wr = 490 fb, R = 11.8, er to t = 38.3 



TABLE II: Efficiencies ex for the heavy neutrino signals with 
(Miv B ,mjv) = (2.5, 1) TeV in the nine samples corresponding 
to £, £' — e, fi,T (rows) and the twelve different same-sign 
and opposite-sign dilepton final states X (columns) separated 
according to the 'low' and 'high' E™ 1BB selection criterium (see 
text). 



Same sign Opposite sign 

Low E% iBB High E^ iBB Low E^ iBB High E™' 1 



££' \ X 


ee 


efi 


/'/' 


ee 


efi 


/'•/' 


ee 


efi 


/'/'• 


ee 


efi 


fifj 


ee 


20.9 








0.4 








21.5 








0.4 








efi 





16.5 








1.8 








15.8 








2.8 





er 


1.6 


1.4 





1.6 


1.3 





1.5 


1.3 





1.9 


1.6 





fie 





11.0 








9.4 








10.8 








9.3 





fifi 








9.0 








8.5 








8.5 








8.8 


fXT 





1.0 


0.9 





1.9 


1.7 





1.0 


0.9 





2.1 


1.8 


re 


0.5 


0.5 





3.5 


3.2 





0.8 


0.7 





3.7 3.3 





Tfl 





0.4 


0.5 





2.9 


2.8 





0.7 


0.6 





3.3 


3.0 


TT 


0.1 


0.2 


0.1 


0.5 


0.8 


0.4 


0.2 


0.2 


0.1 


0.5 


0.9 


0.4 



fb; (2) M Wr = 2 TeV, m N = 1 TeV, for which <j Wr = 
1580 fb, R = 14.2, <7tot = 104 fb. For each case we 
generate nine signal samples with £,£' = e,/x, r. These 
samples allow us to obtain the signal resulting from a 
heavy neutrino with arbitrary mixing, by an appropriate 
averaging with weights ww — \vf N \ 2 \v^ N \ 2 . (Note that 

J2 w w = 1 because l w «vl 2 = !•) The efficiencies Sx 
after cuts for the nine samples are given in Tables HJ and 
IIIII for the mass points 1 and 2, respectively. Notice that 
the Ex are rather similar for the two cases, with some 
migration between the low E™ lss and high i5™ lss final 
states due to the different kinematics. Then, the number 
of signal events for a given final state X is 

Sx = fftotL ^2 w w £ x ■ ( 16 ) 
w 

with er t ot given by Eq. ^ . 

We investigate the prospects for measuring the heavy 
neutrino couplings by simulating pseudo-experiments 
with selected mixing benchmark points given in Table lTVl 
and assuming that the number of observed events in the 
pseudo-experiment equals the signal Sx (which differs 
from one mixing benchmark point to other) plus the SM 
background Bx- Then, we extract the signal cross sec- 
tion and the heavy neutrino mixing parameters from the 
pseudo-data by minimizing the quantity 

2 _ ^ (AgP - Nff 
x — M ex p ' ( ' 

where 

V^ xp = S x + B X , (18) 
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TABLE III: The same as in Table E2 for (M WR ,m N ) = (2, 1) 
TeV. 



Same sign Opposite sign 

Low E'£ iss High E'£ iss Low E!£ iaa High Eg lss 



U! \ X 


ee 


f:'/i 


fifi 


ee 


f:'/i 


fifi 


ee 


eft 


mi 


ee 


e/i 


fifi 


ee 


21.8 








0.2 








21.8 








0.3 








efi 





18.1 








1.3 








17.3 








1.7 





er 


1.8 


1.6 





1.3 


1.1 





1.7 


1.5 
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1.4 





fie 





15.3 








5.4 








14.9 








5.4 





/'/' 








12.7 
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12.4 
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re 
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2.6 
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2.8 
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TT 
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FIG. 1: Limits on unitary heavy neutrino mixings (k = 1) at 
68.3 % CL for the two mass benchmark points (Mw R , mtf) = 
(2.5, 1) TeV (light green) and (Mw R ,m N ) = (2, 1) TeV (dark 
green) in the (v^ N ,vj^ N ) plane (top) and in the (v^ N ,v^ N ) 
plane (bottom). In the top panel, the dashed lines delimit 
the 68.3 % CL regions obtained without separation of the 
dilepton samples into 'high' and 'low' E^' ss . 



ues used in the pseudo-experiments (see Table IIV|) and 
the solid regions are the 68.3% confidence level (CL) 
limits on the two independent normalized mixings v^ N 
and vff N extracted from the pseudo-experiments, as dis- 
cussed above. Notice that the normalized mixings satisfy 
\ v eN\ 2 + \ v ^n\ 2 + \ v ?n\ 2 = 1 by definition, so the region 
outside the unitary red arc is not allowed. The dashed 
lines correspond to the 68.3% CL limits that would be ob- 
tained without splitting the dilepton samples into 'low' 
and 'high' E™ lss (for each mixing benchmark point, the 
pair of outer lines corresponds to My/n — 2.5 TeV and 
the inner lines to Mw R = 2 TeV). In most cases, the de- 
generacy that we have discussed qualitatively in Eqs. (fT3)) 
is manifest. The only exception occur for points A (G) 
where the absence of electrons (muons) in the signal im- 
plies a mixing with the muon (electron) close to maximal. 



TABLE IV: Mixing benchmark points considered in the 
pseudo-experiments. 



Label 


V R 




V R 


A 





1 





B 


0.3 


0.95 





C 


0.3 


0.7 


0.65 


D 


0.5 


0.5 


0.71 


E 


0.7 


0.3 


0.65 


F 


0.95 


0.3 





G 


1 









is fixed for each benchmark point, and 

Nf (5-tot , €n , tf N ) = otot L ^2 w w [y% , t$v ) eg 

w 

+B X ■ (19) 

The above quantity is a function of the two independent 
flavour mixing parameters v^ N , v^ N and the total cross 
section a tot (we use the hat notation to distinguish the 
fitted values from the 'true' ones used to generate the 
pseudo-data sample). The inclusion of systematic uncer- 
tainties is not essential for our analysis, which is more 
focused on the method to extract heavy neutrino mix- 
ings from data, rather than on the precise value of the 
limits obtained. In any case, the statistics are moder- 
ate for the two mass points considered and, therefore, 
we do not expect a dramatic change of the results pre- 
sented here when systematic uncertainties are properly 
accounted for. 

In Fig. Q] we present our results for the two mass 
benchmark points considered. The upper plot shows the 
limits on the two independent normalized mixings v^ N 
and v^ N . The white dots correspond to the input val- 
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FIG. 2: Limits on the heavy neutrino mixings at 
68.3 % CL for the two mass benchmark points (Mw R , ttin) = 
(2.5, 1) TeV (light green) and (Mw R ,m N ) = (2, 1) TeV (dark 
green) in the (V^, V^) plane (top) and in the {V^ N ,V^ N ) 
plane (middle). In the lower panel, the limits on the normal- 
ization parameter k of Eq. Q are shown. 



presented before for each mass point. The absolute mix- 
ings Vin equal the normalized ones (in Fig. [TJ multiplied 
by ft, which is obtained from the fitted cross section a toi 
and Eq. (fl"2")) . As in the previous case, the allowed regions 



for the mixings at 68.3% CL are shown in the (Vjjj, V^ N ) 



and {V^ N: V^ N ) planes (top and middle plots, respec- 
tively). The corresponding limits for the parameter k 
are given in the lower panel (color codes are the same as 
in Fig. [T]). We observe that, provided the statistics are 
good and the systematic uncertainties for the total cross 
sections (which we have ignored) are controlled at the 
~ 10% level or better, the unitarity of the mixings can 
also be tested. 



IV. DISCUSSION 

In case new physics is discovered at the LHC, it will be 
essential to determine its properties. Here, we have ad- 
dressed one particular example that represents a variety 
of new physics models related to neutrino mass genera- 
tion: the production of heavy neutrinos in the context of 
a left-right symmetric model. We have investigated the 
determination of the couplings of a heavy neutrino N 
to the charged leptons, showing that an inclusive event 
counting is not sufficient for this purpose. Basically, this 
happens because the electrons and muons in the final 
state can be produced promptly in the N decays, or from 
secondary decays of r leptons produced in the decays 
of the heavy neutrinos. We have pointed out that one 
possible solution for this degeneracy is to discriminate 
the dilepton events based on the missing energy, so that 
prompt electrons and muons can be distinguished from 
the secondary ones. 

Obviously, this kind of degeneracy is not a specific fea- 
ture of this model, but a generic one in any SM extension 
with new particles coupling to the charged leptons. The 
most immediate example is type-I seesaw without left- 
right symmetry, in 



pp -> W* -> £N -> Wjj , 



(20) 



although in this case observing the signal is very de- 
manding let alone measuring the heavy neutrino 
couplings. More interesting is the production of heavy 
neutrino pairs via a new Z' boson [15l4l8j . 



PP 



Z' NN -> ii'jjjj . 



(21) 



Similarly, in type-Ill seesaw one has dilepton signals re- 
sulting from the production of the new heavy leptons 
E,N 



The lower plot presents the same benchmark points and 
resulting limits but viewed in the (v^ N , vj} N ) plane. 

In Fig. [5] we repeat the analysis for the same bench- 
mark points assuming that the measured value of <Jw R 
corresponds to gn — and equals the theoretical value 



pp -> EN -> Wjj j j . 



(22) 



In all these cases, the procedure presented here can be 
straightforwardly applied to determine the couplings of 
the heavy leptons. In case more than one heavy neu- 
trino is produced at the LHC, the method presented here 
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can be generalized by introducing appropriate kinemati- 
cal cuts to single out the contribution of each heavy neu- 
trino. The practical implementation of this sophistica- 
tion (which would be necessary if several heavy neutri- 
nos were discovered) strongly depends on the particular 
values of the heavy neutrino masses, and falls beyond the 
scope of this work. 
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